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Abstract 

We present scattering from many body systems in a new light. In place of the usual 
van Hove treatment, (applicable to a wide range of scattering processes using both 
photons and massive particles) based on plane waves, we calculate the scattering 
amplitude as a space-time integral over the scattering sample for an incident wave 
characterized by its correlation function which results from the shaping of the wave 
field by the apparatus. Instrument resolution effects - seen as due to the loss of 
correlation caused by the path differences in the different arms of the instrument 
are automatically included and analytic forms of the resolution function for different 
instruments are obtained. Each element of the apparatus is associated with a corre- 
lation length (or time). These correlation lengths, determined by the dimensions of 
the apparatus are generally much smaller than these dimensions and larger than the 
wavelength. As is well known, these are the conditions for the validity of geometri- 
cal optics so that the conventional treatment, where the scattering is calculated by 
the van Hove plane wave approach and the trajectories through the instrument are 
treated classically, is usually valid. In the present approach analytic expressions for 
the correlation functions are obtained. The intersection of the moving correlation 
volumes (those regions where the correlation functions are significant) associated 
with the different elements of the apparatus determines the maximum correlation 
lengths (times) that can be observed in a sample, and hence, the momentum (en- 
ergy) resolution of the measurement. This geometrical picture of moving correlation 
volumes derived by our technique shows how the interaction of the scatterer with 
the wave field shaped by the apparatus proceeds in space and time. Matching of the 
correlation volumes so as to maximize the intersection region yields a transparent, 
graphical method of instrument design. PACS: 03.65.Nk, 3.80 +r, 03.75, 61.12.B 



1 Introduction 

The usual treatment of scattering from many particle systems assumes the initial 
beam to be represented by a plane wave and shows that the scattering cross section 
is proportional to the Fourier transform of the correlation function of the density 
fluctuations. For inelastic (energy changing) scattering this is the van Hove time 
dependent correlation function, in the case of elastic scattering it is the static pair 
distribution function |l|. The usual treatment then integrates the cross section for 
monochromatic plane waves over the spread of momentum in the incident and scat- 
tered beams. 

In this work we show that another approach is both possible and interesting. As 
we do the entire calculation in real space-time our approach sheds new light on the 
emergence of the Fourier transform in the scattering cross section. A crucial feature 
of our approach is that the influence of beam preparation on the resolution of the 
measurement is expressed by the correlation volume of the incident beam, which in 
turn, is determined by the optical (in the sense of normal 12] and time dependent M 
optics) properties of the beam preparation elements, e.g. slits and choppers. 

We will see that the lateral width of the beam correlation volume determines 
the possible spatial extent of the correlations which can be observed in a given ex- 
periment, and is inversely proportional to the q resolution, as in the conventional 
treatment. The longitudinal extent of the beam correlation volume or the correla- 
tion time of the incident beam likewise sets the limit to the system correlation times 
which can be observed. The correlation volume Vc is the conjugate volume to the 
momentum space volume Vp and VpVc « h!^. 

In a previous work y we have discussed the neutron spin echo spectrometer 
[g from a similar point of view. We have shown that in that case the equivalent 
correlation volume is split longitudinally by the 'spin echo time', tnse, which can be 
much longer than the correlation time of the incident beam. This splitting of the 
correlation volume results in the fact that the spin echo spectrometer measures the 
time dependence of the scattering system correlation function directly. When the 
incoming beam correlation volume is continuous (not split), then the scattered wave 
is an integral over all space-time separations falling within the correlation volume of 
the incident beam and the result is the Fourier transform of the space-time correlation 
function of the scattering system as will be seen below. 

1.1 Correlation functions and scattering 

The concept of correlation (coherence) volume has been applied to visible light optics 
for more than 50 years including such problems as the effect of the illumination on 
the resolution of a microscope . The correlation properties of a beam as a result of 
the optical properties of its preparing devices are known as the van Cittert-Zernike 
theorem [0, and |7| and we will see how this can be applied to both time - 
dependent and -independent scattering. 

A frequently used technique for deconvoluting the instrument resolution from the 
measured data is to use the Fourier transform of the resolution function B . In the 
present work we show the physical significance of that Fourier transform and show 
how it is determined by the optical properties of the beam. It has been shown in |p| 
that this Fourier transform may contain more information than e.g. the width of the 
resolution function. In fact it is the correlation function at long distances and times 
that is often of the greatest importance. An additional feature of our treatment is 
that the instrumental resolution is seen as an integral part of the experiment not 
something that is put in ad hoc after the calculation. 
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It is well known that scattering experiments depend only on the spectrum 

of the incident beam and hence can be considered as the incoherent superposition 
of the results of experiments done with plane waves distributed according to this 
spectrum. This is easily seen in the context of the present work where we show 
that the scattering cross section depends on the auto-correlation function of the 
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incident beam. The Fourier transform of this auto-correlation is 
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In ad- 



dition it is usually assumed that sources produce 'chaotic' beams, i.e. beams with 
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A ( fci j A* ( ^2 ) ) = CJ f fci - ^2 ) , d] , so that in any case it is only A( kj 

which is significant. This is the case for 'stationary' beams where ( \ip {x)\ ) is inde- 
pendent of X. 

These points, and the possibilities of measuring (yl ( fci ) A* ( fc2 ) ) have been 



discussed in 1 10 and the references therein, and have been emphasized by Mezei ||ll[ . 
For these reasons we feel that it is more convenient to talk about correlation volumes 
rather than coherence volumes as has been the practice in optics [^ since it seems 

reasonable to reserve the term coherence for those cases where ( ^ ( fci 1 yl* ( fc2 1 ) 
plays an important role. 

The correlation length Sc discussed here and in section has real physical mean- 
ing (as of course does its equivalent the spectral width Afc). This can be seen by 
considering that in order to observe the famous two slit interference pattern (Young) 
it is necessary to illuminate the two slits with correlated radiation. This is normally 
produced by having the two slits located in the main peak of the diffraction pattern 
produced by a third slit up-stream from the pair producing the interference pattern. 
The van Cittert- Zernike theorem is the rigorous expression of this point. The two 
point correlation function can be measured with such a double slit experiment by 
measuring the modulation (visibility) of the interference pattern as a function of the 
slit separation. Of course the physical situation can be described according to the 
spectral viewpoint, where the incident beam is considered as an incoherent superpo- 
sition of different k vectors, by considering that for increasing slit separations the 
interference patterns for the different incident k vectors begin to cancel out. 

Recent experiments involving Mossbauer scattering of synchrotron radiation have 
been described in terms of the beam correlation volume p2[ . These experiments can 



also be interpreted in terms of 



A k 



However the authors have shown that 



the concept of a coherence (correlation) length is very useful for understanding this 
exciting experiment. 

In the following we shall see that these ideas have a rigorous basis. We have 
presented the main features of these ideas with some physical examples in |0] . 

1.2 Examples of correlation lengths and times 

We now discuss two simple examples to show the magnitudes of the correlation 
lengths (times) introduced by the shaping of the wave fields produced by the optical 
elements of the apparatus. As formulated by Marathay "..the field from a non- 
coherent (uncorrelated) source acquires coherence {a finite correlation length) by the 



very process of propagation." (see |14 , p. 105.) This is nothing else than what occurs 
when a random 'white noise' signal acquires a finite correlation time as a result of 
passing through a low pass filter, something which can easily be observed with an 
oscilloscope. 

For elastic scattering, a beam collimated with slits of width a separated by a 



distance L, will have an angular width of A^ ~ a/L and a spread in momentum per- 
pendicular to the beam of A/cj, ~ k^AO ~ k^a/ L. This will result in an uncertainty 
in momentum transfer q of Aq ~ Aky. 

On the other hand, seen according to quantum mechanics, a beam going through 
a slit of width a will have an uncertainty in momentum Ak' r^ 1/ a and this will 
result in a lateral spreading over a distance 5 ^ Ak' L/k^ — L/kxa at a distance L in 
the X direction. We see that the correlation length i5 ~ 1/Ag so that both arguments 
will lead to the same resolution. In order to put the argument on a firmer footing we 
need to show that a wave acquires a correlation length 5 on passing through a slit of 
width a. 5 is the length associated with a Fraunhofer diffraction pattern of a slit. The 
problem is then to show that this length can arise in the case of Fresnel diffraction 
(finite distance L). This is the content of the van-Cittert Zernike theorem (g, |^, 
0) mentioned above. In this way we can describe the scattering process in terms of 
diffraction at the defining slits and the correlation length of the incident beam. We 
will be able to describe the scattering in real space instead of in the Fourier transform 
space defined by if . 

For a small angle scattering instrument with wavelength of 10 Angst, and coUi- 
mating slits 2cm wide separated by L = lOm the correlation width will be about 
1600 Angst., i.e. much larger than the wavelength but much smaller than the slit 
width. Correlations over larger distances than this cannot be observed with the 
assumed instrument. 

Applying the same argument to a chopper for massive particles we see that 
for a chopper opening T and time of flight to we will have AE/E = 2T/to = 
Alo/lo, {hw — E = mw^/2) by the usual classical arguments. Quantum mechani- 
cally, according to the uncertainty principle, passing through a chopper of width 
T will result in a spread of energy given by Auj' ~ 1/T and this spread in energy 
will result in a spread in arrival times At/to ~ Av /v ~ Auj' /2uj ~ 1/2loT . Thus we 
should expect a correlation time of At ~ to/2LoT ^ 1/ Auo. For a time of flight instru- 
ment with pulse width T = 10~^s and a flight path of L = 5to, neutron wavelength 
A = 5 Angst, we obtain At ~ 30ps and this is the limit to the correlation times which 
can be observed with such an instrument. 

2 A new look at scattering from many particle systems 

We will now show how the effects discussed above are contained in a quantum- 
mechanical calculation of the scattering. Following the usual treatment of scattering 
P, sec. 2.1, |l3] the scattered wave is given by 



i)sc (T", t)^ j Go{l^ ~rt,t~ ts) V {rt, ts) Ac [rt, t^) (fr^dts (1) 

where ipse (j^) is the exact solution for the wave function at the position rt, and Go is 
the Green's function for the unperturbed problem (V = 0). The notation is explained 
in figs. 4 and 6. As is usual in van Hove scattering we make the Born approximation 
by replacing ipse (rt) by the incident wave at the position of the scattering sample, 
An i^) and take (for the case of neutron scattering) 

V r?, t = V hS {rt-rtit)) = bp r?, t) 2 

^ — ' m m 

i 

(for coherent scattering), i.e. the Fermi pseudo potential where Tf is the position 
of the «*'* nucleus, the sum is over all nuclei in the scattering system and p {r^, t) as 



defined in (0) is the particle density of the scatterer.R 



^- (^' t) = jGA^-rt,t-u)v in, t.) A. irt, Q d^rA (3) 

The integration is carried out over the sample volume and the time interval during 
which the sample is exposed to the beam. 

Equation (0) describes the scattered wave as a superposition of waves produced by 
many scattering events each occurring at the different space time points (r^ , ts){Fig. 
1). This is in contrast to the usual van Hove treatment [n], where one begins with 
the Fermi golden rule for time dependent perturbations and shows, replacing the 
energy conservation delta function by its Fourier transform, that the probability of 
scattering is given by the space-time Fourier transform of the correlation function of 
V {rt,ts). In the appendix to jj] we have shown how one can derive the van Hove 
result in a more direct way. Introducing rt — rf = 6 , ts — t'^ = t we have from (g|) 
for the intensity at the detector, {rd,td)'. 



|V'sc(r?,td)|^= d\sdts I d'^r'Jt'^Goiri -rt,td~ts)Gl{ri -7lM-t')j x 

X V (r?, U) V* (;t, t;) V™ [rt, U) i,*^ (^, t^) 

= / d^Tsdts / d^SdrGo {r^ - r^, trf - t,) Gl (t^ - r^ + T ,td ~ ts + t) 

X V {rt, ts) V* (rt ~ 5 ,ts - t\ ijjin {rt, t^) V'*„ (jt - S ,ts ~ t\ 

(4) 
Substituting the Fermi potential we find 

\^scijt,td)\^ = j d^Tsdts / d^SdrGo {ra - rt,td - ts) G* (rd - rt + S ,td-ts+TJ 

(5) 
xGs(l',T)i?„ (l',r) (6) 

Here 

Gs {T, r) = (^p in, ts) p(n-T, ts - r) )^ (7) 

is the (van Hove) density-density correlation function and is independent of rj, ts for 
the usual case of a homogeneous system. The result depends on the pair correlation 
function iQ) due to the fact that the intensity is proportional to the square of the 
amplitude as can readily be seen in the above derivation. Because we restrict our- 
selves to the Born approximation the result is independent of higher order correlation 
functions. 



Ri, 



(^, TJ = {i^in {rt, ts) 1p*n [rt -'S^,ts-Tjy (8) 



^Iii the interests of simplicity we omit prefactors and constants in the body of this work con- 
centrating on those terms which are essential to the physics. We give a complete treatment in the 
appendix where we also show how to extract the cross section from our results. 



is the auto-correlation function of the incident beam (See fig. 2 for the time indepen- 
dent case.). The brackets 0^, and ()j„ indicate statistical averages over the sample 
and incoming beam ensembles respectively. 

Rin ( '^ I ■'' ) 'will bs seen to be responsible for the contributions of the incident 
beam to the resolution of the measurement, as, according to (0), the intensity only 
depends on the correlation function Gs at values of f (5 , r 1 where Rin is significant. 

The wave function incident on the sample V'm, is given by 



i^in (rt, ts) = / (frodtoGo {r^ -rt,ts- to) tpo (r^, to) (9) 

where ipo {r^ , to) is the incident wave function at the entrance of the apparatu^ We 
then find 



" (^' ^) " / '^''^°^^°^° (^' ^« " *°) ^° (^ -T,h- r) (10) 



Ri 

where rf = rt —r^ ; ti = ts—to and we have made the assumption that the incident 
beam is completely uncorrelated, i.e. {ipo {r^, to) V'o (^'j *o))m = ^^ (^ ~ ^') ^ (*o — t'o) 
as is usual in the derivation of the van Cittert-Zernike theorem H , jq] and M . By the 



Wiener-Khintchine theorem H we see this is equivalent to 
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(t) 



const, i.e. 



the incident spectrum is broad compared to that selected by the instrument. 

In a scattering experiment we measure the integral of \ijjsc{rd,td)\ over the 
entrance area of the detector and, in the case of a time of flight experiment, a definite 
time interval, dtd. From (o) the intensity at the detector, I^, is given by 



d^rddtd\ipsc{rd,td)f ^ d^Vsdtg d^SdrRout (S , r, fj ,<2J G^ f (5 ,TJi?.„if(5 ,t, rf,iij 

(11) 
where 

Rout (T, T, rj, ta) = / d\ddtdGo (rj, ts) Gl (fj + T,t2 + r) (12) 

and r2 ^ rd - rt , t2 ^ td - ts. 

Although ( [ll| ) contains an explicit dependence on only space and time variables, 
the dependence on the energy and momentum transfer will be seen to emerge from 
the variation of the phases of the Go functions. 

Now we see there is a symmetry between the integral over the detector parame- 
ters, {vd, td), of the product of Green's functions depending on (rj , t2). Rout [5 , r j , 
(equ. |l2| ) and the integral over the input slit (r^, to) of the Green's function product 

^Strictly speaking we have ipini^ ,is) = J cProdto — ° a-^' ° — —i)o(rS,to) where Gois a 
Green's function satisfying the boundary condition Gn = on the surface of integration, and d/dTT' 
is the gradient in the direction of the surface normal [t8|,|3|. However this only introduces a change 
in the prefactor. See the appendix for a complete treatment. 



depending on {ri,ti), Rin \S ,r), (equ. |lO| ) so that we only have to do the cal- 
culation of one of the pairs, say the latter. We then obtain the former integral by 
applying the symmetry rules: 

5 -^ -6 , T ^ -T, t2 ^ti, r2 ^ ri, rt ^ -rt, tg -^ -ts. (13) 

While functions such as Rin,out ( ^ i''") ^^^ usually called 'correlation functions' 
they really represent a loss of correlation as their arguments increase. In the present 
context it is instructive to keep this in mind by referring to them as '(loss of) corre- 
lation functions'. 

Thus while i?i„ I 5 , t, r]*, ii ) represents the loss of correlation between pairs of 

rays starting at the entrance slit ("P'o), propagating to neighboring space time points 
in the sample (rj,ts), {rt',t'g), Rout (^ i''') represents the loss of correlation due 
to the difference in path lengths in propagating from two neighboring points in the 
sample (rjjis), {^f ,t'^) to the detector {T*d,td)- GsiS ,rj is the loss of correlation 

between the points (r^,ts), {r^',t'^) due to the internal dynamics of the scattering 
system (See fig. 3). Implicit in the above discussion is the fact that the detector 
responds to the value of \ipsc\ at a single point integrated over the detector area. 

The above discussion brings out the fact that the scattering process is intrinsi- 
cally an interference phenomenon, the cross section being a measure of the loss of 
correlation between waves following paths separated by I 5 , r j at the sample. 

In this paper we limit ourselves to pure coherent scattering. The separation of 
the scattering into a coherent and an incoherent part can be carried out exactly as in 
the conventional treatment. ||], |^. Also in the interests of simplicity we choose to 
work with a real 'hydrodynamic' density function p(T*,t). Replacement of this by 
Si^C/^ — T*i(i)) yields results consistent with the conventional approach (l), (see 
also [|). 

3 Elastic Scattering in Real Space 

In this section we specialize to the case of scattering from a static system. We will 
consider time dependent systems in the following sections. 

Using the known form of the Green's function for the time independent Schrodinger 
equation, 

we calculate the input product of Green's functions 



by expanding 



r s - r „ - r , - r „ = 



ik{\T*^-ri\-\T*'^-r^\) r-^^-. 



S ■ (T^s - T^o) S ■ {T*s - ~r*o) S ■ {T*s - ^o) 



Ll\ 



(15) 



T*i — T*s is the location of the center 
of the sample and we rewrite ~r^s =^ X^s + ~^s where ~e^s is now a small quantity 
(see fig. 4). Then 



Rin I 1 = e --i I d Tot --i = e --i / rfyoC "i = 

where we have taken a one dimensional sht of width 2a and k i = k T*i/ri and 
r^ = |T*i|. (5y is the component of 6 perpendicular to the incident beam (see fig. 4). 

The result ( |l6|) is known as the van Cittert-Zernike theorem, f^li O, and 
has the form of the Fraunhofer formula. Note that this formula is applicable to 
geometries where one would not normally observe a Fraunhofer diffraction pattern, 
except in the case of a spherical wave converging on the point of observation. In 
applications to optics the phase factors in (|l6|), depending on the position of the pair 
of points with respect to the center of the sample, as well as on the distance between 
them are often neglected so that the correlation functions can be considered as real 
functions depending only on the distance between the two points. However the phase 
factors play a crucial role in scattering, bringing in the influence of the sample size on 
the resolution and the dependence on the energy and momentum transfer, as shown 
below. 

By symmetry ( [l3[) we can write 

Rout{T)=e-^i^^'e"'^^2di^jr^\ (17) 

where 6' represents the component of 6 in the direction perpendicular to the direc- 
tion of the scattered beam, fc 2 = kT^2/r2 and d is the width of the detector slit (see 
fig. 4). 

Then writing the steady state form of (O) we have 



Id^ f d'cs f dHRout (^, n) Gs (t) R„, {t. n) (18) 

Substituting (|6|) and (|l^) in (^ we have 



d\s I dHe-'^-^e'''^^2d \ ""^r ? \gAT] 2a 



smA:^— \ /_.\ / smk^— \ 111 



j jv"^ / n / \ fc^ / ' r' ri r2 

r2 ' \ ''1 / 

(19) 

The term exp ( ik '^' J comes from the phase factor which is normally neglected 

in optics. {Note that if — k2—ki,a point which will be discussed in the next 
section.) In the present case this factor means that moving the location of the 
correlation regions defined by the two svax/x functions (varying ~t s) in ( |l9| ) causes 
an additional difference in the two optical paths going from the entrance slit to 
the detector slit via the two points separated by (5 , i.e. the two paths in fig. 4. 
This optical path difference represents an additional loss in correlation as the pair of 
points separated by a constant 5 move away from the optical axis and determines 
the influence of the sample slit width (26) on the resolution. 



3.1 Small angle scattering 

For simplicity we will now consider the case of small angle scattering where the only 
significant component of 6 is Sy = ys — y'^ and 6y « Syi . 

Performing the Cs integration in ( |l9|) over a slit of width 26 ( J_j^ de^ j defining 

the sample size yields 



Id^ I dHe-'^-^Gs 



6 \2d 



sinfc^^ 

r2 

, (5„,d 



26 


[sinfeSM 


2a 


sin/c^ 



(20) 



where we have neglected the extension of the sample in the z direction. 

For a, 6, d small, corresponding to a high 'if resolution, the functions in square 
brackets approach unity for the values of 5 for which Gs is significant, and we have 
the usual result 



5(T) 



f d^Se-'^-^Gs(t) 



(21) 



The same result is obtained if we take the incident and outgoing states as plane 
waves. The result of an experiment for non-negligible a, 6, d must then be found 
by taking the convolution of ( |2l| ) with an instrumental resolution function. In the 
approach of ( pO| ) however, the resolution is contained in the (loss of) correlation 
functions in the square brackets. Equ. ( po| ) clearly shows how correlations between 
points in the sample separated by larger values of S contribute less to the intensity. 
However due to the nature of the sin x/x functions there are some contributions to the 
scattering from these larger values, and some regions of S contribute with a negative 
sign. These effects, at the 10% level, which do not appear if we approximate the 
sin x/x functions by Gaussians of an appropriate width as is often done in practice, 
become more important for measurements at poor resolution such as in quasi-elastic 
line broadening. 

3.1.1 Comparison virith conventional treatment 

We see that equ. ( po[) is the Fourier transform with respect to S of a product of 
G I S ] with a function H (a, 6, 6y) which is itself the product of three functions, /i2 {d, 5' ), 
h' (6, Sy)a.nd hi (a, Sy). Thus the result ( po| ) is the convolution of the Fourier trans- 
form of G ( 6 j, i.e. S (if), with the Fourier transform of H. The latter is the 
convolution of the Fourier transforms of each of the h functions 
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(22) 



(ri/fc is a normalizing factor so that rja (qy) is normalized to 2a and u (x) is the unit 
step function). In rii,{qy) the term (a/ri) is replaced by {b/r') and in r]ii{qy) it is 
replaced by {d/r2). The convolution of two rectangular functions in qy, rja and rj},, 
yields a trapezoid with a base bounded by 



\qy\<kia/n+bil/n + l/r2)) 



and a top bounded by 



\qy\ <fc|a/ri-6(l/ri + l/r2)|. 



just what is expected on geometrical grounds for a collimator consisting of two slits 
of width 2a and 26 separated by a distance x as seen by a point r^ away from the 
second slit (see fig. 5). The subsequent convolution with rjd{qy) yields the effect of 
finite detector width. From (Efl) the correlation length of the input function hi is 
seen to be j/c = ri/ka, while (E2h gives the width i/q in q, to be j/g — ka/ri. Hence 
He ■ Vq — ^ or He ■ ijp = h, where the momentum width j/p — hyq. 

Thus our result (PQ ) is identical with the usual treatment which yields the convo- 
lution of S i'q*) with the resolution function but our view of the scattering is quite 
different. We see the scattering process as the interaction with the sample of a state 
whose wave function has a certain spatial auto-correlation function determined by 
the optical properties of the devices used to define the incoming beam. This has the 
effect that only pairs of scattering points whose separation is less than the correlation 
length of the incident beam contribute significantly to the scattered wave and pairs 
with different separations S are weighted by the auto-correlation function of the 
incoming beam evaluated at S . Analogous statements apply to the outgoing beam. 

Similar considerations have been used by many people {e.g. [g|) as a means of 
deconvoluting the instrument resolution from the measured data. However, in the 
work cited no mention is made of the physical significance of the Fourier transform 
of the instrument resolution function (our H (S) is called R{k) in (Q)) as the (loss 
of) correlation functions of the incident and scattered wave functions. In the above 
treatment we see how the resolution function is the result of the optical properties 
of the devices defining the incident and scattered beams. 

3.1.2 Discussion of the result for small angle scattering 

In the previous section we have seen that the contribution to ( jV'sc (^)l ) from the 
correlation function of the incident beam represents the effect on the resolution due 
to the coUimation of the incident beam and the correlation function of the scattered 
beam results in that part of the resolution due to the angular resolution of the 
scattered beam. Each element of the optical system defining the beam is seen to 
contribute a factor to the overall correlation function. 

In the following we try to illustrate the relation between our treatment of the scat- 
tering problem and the conventional approach by concentrating on the contribution 
of the incident beam to the resolution. If we write 



*mc(r?)- / d-'kiAl ki)e'''-^' (23) 



then we have 



(* „, (r?) n^c (^) )^ = / d\s f d^h f d% I^A [ki) A* (k^) ) e^^ 



■rt—k-, -r^ 



'^U \ A (1. \\^ 



ih\ ■ I rt — r^ 



d'ki{^\Aiki)\')e V V (24) 

Substituting this into ( p^ ) and writing ipout — e^^^'"^ we obtain 

J d\s J d^h(^\A{ki)\^) f d'Se-'^-^G, (t') ^ Jd^ki(^\A{ki)\^)s{lt) (25) 

In this case we have taken if — k2 — ki and we obtain S (if) averaged over the 
momentum spectrum of the incident beam. 



In the previous section we took if — fc 2 ~ ^i , *• e. the difference between the 
nominal final wave vector, and the nominal incoming wave vector. By this we mean 
the wave vector corresponding to a central ray of the incident or scattered beam. 
Then, in that case, the distribution in fci, or what contains the same information, 
the correlation function of the incident beam produced by the optical properties of 
the defining slits, gives the usual result of S (if) convoluted with an instrumental 
resolution function. On the other hand, in this section by taking 



-t=k2- ki (26) 

we have assumed perfect knowledge of fci , fc2 ^^nd we complete the calculation (equ. 
E9) by integrating over all values of if corresponding to the spread in ki values in the 
incident beam. This will of course yield the same results as the previous treatment. 

3.2 Elastic scattering at arbitrary angle 

We will assume a cylindrical shaped sample (radius Rg) with its axis perpendicular to 
the scattering plane and thus consider the problem as a two dimensional one. Then 
the integral over the sample volume in (pjj[) can be written as 



d^Ese'^—^ = I I esdes / dOe'''^^^^ = I / e,de^27rJo ( k^ 



JO JO 

2' 



2^^s l,_l ' l (27) 

r' 

where / represents the dimension of the sample perpendicular to the scattering plane. 
The final result is then 



/' 



d^ra li^sc iri)f = f dHe~'^-^Gs (T) H (T) (28) 



where 



/sinfc^\ /jjfc«k)\ /sinfc^\ , 

is the Fourier transform of the resolution function. The terms depending on a, d 
are identical to those discussed above in connection with small angle scattering (EQ). 
The Fourier transform of the term containing i?s is given by 

JsdddOe^^'-' (^^) - 2nJdSJo {qS) 0^) = (30) 

lo {q>4-} 
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3.2.1 Discussion 

Just as in the small angle case the result in the general case (^8|) is seen as the Fourier 
transform of the product of G^ with a product of functions, each representing the 
contribution to the resolution of an individual element of the beam handling system. 

Corresponding to the sin x/x functions in one dimension the function Ji {x) jx in 
two dimensions represents a square shaped resolution in q but in the latter case the 
resolution is a function of |"g*| rather than a component of ~<f in the one dimensional 
cases. 

The term involving (5j,a) gives the influence of the width of the entrance slit (2a) 
on the resolution; it depends on r\ but not on r2. The term in {6yid) gives the 
influence of the detector slit width {2d) on the resolution and depends on r2. The 
remaining term gives the influence of the sample size and depends on r' because an 
uncertainty in position at the sample effects the angular resolution in both arms. 

The argument of each term can be written as (S/dd) where the S^ are correlation 
lengths corresponding to the size of the correlations that can be observed with a 
given installation. When any of the terms is small, so is the contribution to the total 
integral, so that the only significant contribution comes from values of 6 for which 
each term in ( p9|) is significant. Thus the resolution is determined by the overlap of 
the different correlation volumes and the measurement can be optimized by choosing 
all regions as approximately equal. 

Equ. (|28| ) has been derived for the case of a single incident energy. It is easy to see 
that if we took an incident spectrum constant in a region ±AK that the result will 
be an additional term sinx/a; with x = AK (Si — 62), where 61^2 are the components 
of 5 along the directions of fci_2 respectively. A more detailed discussion including 
the influence of the scattering angle will be presented in a forthcoming paper. 

3.3 Neutron spin echo for elastic scattering. 

In neutron spin echo for elastic scattering, the relative phase shift between the two 
spin states is given by if ■ (5nse with if = k 2 — fci H- Then the average beam 
polarization will be given by equ. (p5|) with S replaced by I (5 — Jnse ) and to this 
order the result is independent of the momentum distribution of the incident beam. 
The dependence of the spin echo signal on the incident beam comes from the fact 
that (5nse is a function of fci . We now show this in more detail. 

In [Bj we have described various forms of spin echo spectrometer from a space- 
time point of view. According to this picture the flrst arm of the spectrometer 
causes the beam incident on the scattering sample to be split in space and time, and 
this splitting results in the scattered beam sampling the scattering system at points 
separated by this splitting. In this section we show how the current view point yields 
the same results. 

We note that we are considering only scattering with if — iqy. In the case of 
spin echo the wave function: 



(7^)= / d^kiA{ki)e'^'-^ (32) 



represents the beam arriving at the sample in the absence of the NSE magnetic fleld. 
Turning on the flrst spin echo fleld (length L) introduces a phase shift which 'splits' 
the beam on the sample Q . Calling the components of the split beam on the sample 
tp2 {rt) we have: 



11 



= vl/,„^(r^±5j,SE/2)e±*'^^^/''° (33) 

Since the measured quantity in spin echo [a^) depends on the cross correlation be- 



tween ij)'^ we calculate this cross correlation for the wave functions (33) 



(^2+ (r?) i^r {pj ) = (* ^nc {ri + <5nse/2) vf *„, i^, - <5nse/2) \ 

= Rin^in {Ts -~r'^+ ^NSe) = Rin-in {5 + (5nSe) (34) 

where Rin^in (S) is the auto-correlation function of the incident beam in the absence 
of NSE field. 

Now following (0) the wave function of neutrons scattered by the sample that 
would reach the detector in the absence of a field in the second coil is: 



^tciri) = / d^r,e-^''^-^p{rt)^f (r?) (35) 



where 



fci 



ko (elastic scattering) and k2 is directed along a line from the 



center of the sample to a point on the detector, r^, so that fco^'d = ^2- ^- Thus 
integrating over d'^k2 will be equivalent to integrating over the detector area. 

When the second NSE coil is turned on, the beam reaching the detector will have 
an additional phase shift fe^^'^^^/'"°e'^^''^^^^^/'^) so that at the detector we will have: 



i^U^) = e'"''^ J dV,e-*'=-^eT''=^^"-/2p(r?)vI/,„,(r?±JNSE/2) (36) 

Then 

(O = (Vd(r^)k.|^,(r?)) =y'd%y■dV,y■dMe^'^■(^"'')e-'^■=*--••• 
• ■ ■ /o ( r?) * ^nc {r: + (5nSe/2) p ( K^) ** ^„c (j^ " ^NSe/2) 

The integration over (fik2 gives S ( r'^ — r^ — (5nse ) so that we have finally 

(a,)= /dV'(p(r?)p(r? + <5NSE)),(|*™c(r? + '5NSE/2)|')^ (37) 

= G (^nse) i?i„-i„ (0) = G ((5nse) 

Thus we see that in the case considered here with the range of available ^2 large 
enough, and the approximation that (5nse is independent of ^2,1 we get the exact 
correlation function, independent of the properties of the incident beam. Of course 
we know this from the classical treatment. 
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4 Scattering from a time dependent system 

In the present section we will show how the ideas presented in the previous section 
can be applied to scattering from systems fluctuating in time. In scattering from such 
systems the energy of the scattering radiation is changed (inelastic or quasi-elastic 
scattering) . For simplicity, in this section we will neglect the position dependence of 
the scattering, considering the scattering system to be concentrated at a single point, 
rt and restrict ourselves to massive particle scattering where the energy is defined by 
choppers. I This is equivalent to calculating the scattering integrated over all values 
of If^. We will treat the case of time and position dependent scattering in the next 
section. 

Under these restrictions equ. O) can be written: 



V'sc (rd, td) = / dtsGo [rd -~rt,td- tg) V (r^, tg) tpm (??, *«) (38) 

for the wave function at the detector located at Yd at time td where 



Go {n -^s,td- ts) = e =«(*^-'=) (39) 

is the Green's function (neglecting an unimportant prefactor) for the time dependent 
unperturbed Schrodinger equation. 

The beam will be taken to be incident on the scatterer through a slit (chopper) 
located at the origin that opens during the time interval ~T < to < T. (See fig. 6 
for a definition of the notation used here.) Analogous to the elastic case above, we 
consider the wave function incident on the chopper to be completely uncorrelated 
for different values of to- That is, we assume that the correlation time of the beam 
incident on the chopper is much shorter than the correlation time that will be imposed 
on the beam by the action of the chopper, or, in other words, the beam incident on 
the chopper contains a much broader energy spectrum than will be selected by the 
chopper system. 

The wave function leaving the chopper and arriving at a point rt at time ts is 
givenby0, |I|, § 

ij^„{rt,ts)^ / dtoe^'^<.^^-^o)^o{to) (40) 

J-T 

and 



(V'™(r?,i,)^*„(r?,t;)) = e-'^-e^'^«»^ / di„e-2^-*o/ii (41) 

J-T 

aiT 

where wi = mrl/2ht\ 

This emergence of plane wave-like states from the time dependent Green's func- 
tion has been discussed in ilql. 



^The use of crystals aa_energy defining elements has been briefly discussed according to the 
present viewpoint in ref. |13| and will be discussed in more detail in a later work. 
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In deriving mnj we have taken 

(V'o(io)V':(0) = '5(t„-o 

This (^l|) is the time analogue of the van Cittert-Zernike theorem (^6|) for matter- 
wave optics. 

Following arguments analogous to those in the previous section we find for the 
intensity at the detector integrated over a time channel of width 2Tj_: 

I dU\i>sc(r2M)? = I d^s [ dre^-^e^'^^^GAr)2T'^^2T,^^^ (42) 
J J J aiT a^Td 

-Td -Ts 

where 

2loit 2uj2T 

Qi = ~ , "2 = ^^, ^ = ai + a2 (43) 

El £2 

and ljj = u}2 — i^i is the mean energy transfer of the scattered particles. There is 
a second chopper placed close to the scatterer which opens for tg—Ts < tg < ts+Tg. 
The integral over d£,s is seen to give 2Ts '^'"^ ° with 

....(|.|) ,«, 

The term e*''^= comes from the phase factor which is normally neglected in optics. 
In the present case this factor has the effect that moving the correlation interval (r) 
(varying ^s) in (^ ) causes an additional difference in the two optical paths going 
from the entrance chopper to the detector chopper via the two points separated by 
T, i.e. the paths 1 and 2 in fig. 6). This defines the influence of the sample chopper 
width {2Tg) on the resolution. 
We then have: 

dtd IV'.c {n, td)\^ = 2Tg2T2Td / dre'^^Gs (r) H (45) 

sina2TdSin/^Ts sinaiT 

H= (46) 

a2Td fJ,Ts aiT 

where we have written Gs (t) for (p ( r^", ts) p* {rt, tg — t))^. the time dependent density- 
density correlation function of the scattering system. 

Thus the scattering cross section is found to be proportional to the Fourier trans- 
form of the density-density correlation function of the scattering system multiplied 
by a function of r which is the product of three functions representing the effects of 
the opening times of the three choppers on the resolution of the measurement. In 
this case the scattering is seen as the interaction of the incoming state with the time 
fluctuations of the scattering system. The wave function of the incoming state has 
an auto-correlation function determined by the (time-dependent) optical properties 
of the defining choppers. Only pairs of scattering events separated by times t, less 
than the correlation time of the incoming wave contribute to the scattered wave and 
their contribution is weighted by the beam auto-correlation function evaluated at r. 
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4.1 Discussion 



The situation represented by ( [45| ) is completely analogous to the elastic scattering 
case (p8|). The opening times Ti^2 influence the resolution by an amount depending 
on ti 2 respectively, while Tg contributes to the uncertainty in the travel times in both 
arms and depends on a weighted combination of the two travel times. The argument 
of each term can be written as r/Td where the Td are different correlation times. 
Clearly the measurement will be optimized when all the Td are approximately equal. 

5 A space and time dependent system 

Rewriting equ. (O) we have for the intensity at the detector: 



(frddtd\ilJsc{rd,td)f = / cfrsdts / d^SdrRout 
with Go given by equ. ^. We recall that 



(S ,T,r2,t2]Gs (S ,TJi?i„f(5 ,r, Ff,t] 

(47) 



and we expand 



6 , t'^^ ts 



(48) 



ri -r^ ~ \rt -ro\ - 2 (5 • (r^ - r^) 



(49) 



t's -to U- to {U - toY 



(50) 



where we have neglected terms of the order r^ and S^ . To calculate the product of 
the two Green's functions in (lin) we need to calculate the quantity 



A,; 



(ts ~ to) {f, - to) 



{ts - to) 



\\rt - r^l^ -2~t -{T^ - rt)j x 



t.-t 



_ \rt -rt\ T 2 5 ■ {rt -rt) 

{ts - tof {ts ~ to) 

Using our usual notation 



o {ts — to) 



(51) 
(52) 

(53) 



ts=^t^ + (s=h + 6 



(54) 



rs ^ r^ + es = r i + es 



we have 



(55) 
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\rt~rt\^T( „4' -<o\ 2 (5' • (r^- f^) / L-to 



A...-^ir^^i-^^J+ - -' (^1-^1 (56) 

and finally 



(57) 

|2 ■ 



2J-(X'i) 2|T'irT\ _, /2(5 2T*iT 



Thus the phase of the integrand of Rm (equ. Q) is: 



where 



— Ai„ = -ct;;^T + _k_;^ ■ 6 + {to - 6) Qfi + (rS' - e]f) • /3i (59) 






2^ if ' — ' ft ii 






Then 



J J Piya aiT 

y=-a -T 

(62) 

where we have specialized to the case of a one dimensional input slit going from 
y — —a to y = a. Z(;2> a) is the height of the slit. Note that fc i is defined to be 

parallel to the x axis so that ( A i ) =0. This result (p^ ) is the generalization of 

the van Cittert-Zernike theorem to the case of time dependent matter wave optics. 
Comparison with the pure spatial case mm and the pure time dependent case (El]) 
shows that these are special limits of the general case (p3). 



By our symmetry rules (13) Rout can be written 



J-Ta P2y'd a2Td 

Det 

(63) 

where y' indicates the direction perpendicular to the output arm (T^2 = T^d — T*.s) 
and 
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a2 



2UI2T d ■ k 2 

k.2 1.2 



k2T "ni 6 



t.2 "-1:2 



(64) 



Performing the integral in mc() over the sample volume and scattering time yields 



Ts 






s Jo Jo 



(65) 



where we assumed the sample to be in the shape of a cylinder of radius Rg. In 



equ. 



^ = a2 - «! 



k I k 



2r 



61 to 



(66) 



and 



P2- /?: 



The two dimensional integral gives 



771(5/1 1\ I k I fc2 

— T 



^1 ^2 



(67) 



27r/ ^6,de,Jo(z.6,) = 27ri?2"^^('^^^) 



i^i?s 



The result of evaluating (|47| ) is then 



j d\ddtdHsc{nM)\^ ^ j dHdTe<^^-^-^)Gs (y 



where 



(68) 



(69) 



ijJ = ^2 ~ '^1^ ~9*= ^2^ A; 



(70) 



2Md 



sin/32y'd\ / sina2'?d\ / sin/^T^X / n2 -^1 (^-R. 



Trf 



a2Jd 



T, 



M^s 



27ri?^ 



i^i?. 



5.1 Discussion 

For convenience we rewrite below 



sinoiT^N / smPiya 



aiTo 



(ilyO. 



(71) 



Qfi = 



Ai 



-L.{r-^,r), /3^ = -^(T-^ir 
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a. = -^.(r-V.r), n = ^(r-lf.r) (72) 



M = ("2 - ai) , 



(32- /?: 



Qfi and (3i represent the influence of the input sht and chopper on the overaU 
resolution, fj, and v the influence of the sample size and chopper while a2 and /?2 
represent the efi'ects of the detector slit and time resolution. 

Each of these variables has a term proportional to r and a term proportional to 6. 
Keeping in mind the Fourier transform (p9h we see that the terms proportional to r 
will give the energy resolution and those proportional to 5 will give the q resolution. 
Thus, for example, the r term in ai represents the effects of the first chopper opening 
time, To, on the energy resolution just as in the pure time dependent case, equations 
( p^ , |4^ ) while the S term represents the effect of the first chopper width on the 
q resolution, i.e. the chopper width results in an uncertainty in fc i that gives a 
contribution to the q resolution. Similar remarks hold for each of the terms in ([72[). 



Note that by definition fci , k2 , = but we display these terms in ([72|) to keep the 
symmetry. 

We see that our approach yields directly the mutual influence of the q and lu reso- 
lutions in a compact analytic form. Since each element of the apparatus is associated 
with a definite factor in H, we see directly the infiuence of the individual elements on 
the overall resolution. Modification of these elements can readily be accommodated 
as shown in the next section. 

If we wish we can follow the more usual procedure of replacing each term by 
a Gaussian function with the appropriate width. This then allows an easy way to 
estimate the overall width of the uj, q resolution (given by the convolution of the 
Fourier transform of each term in H) as the square root of the sum of the squares 
of the widths of each term and shows that the optimum is when all the widths are 
equal. 

Each factor in (uW is associated with a certain velocity in I S ,t] space so that 
the relevant correlation volumes can be considered as moving as a function of r. For 
example the volume associated with a2 is moving with a velocity ~w*2, the nominal 
velocity of the scattered wave. Thus we see that the coherence volume associated 
with the detector parameters can be considered as consisting of the region limited by 
■iihL^/md in the 5' direction and ±t2/'^d,k2 in the S^' direction, moving (as function 
of r) with velocity V2 in the x' direction (direction of the scattered beam). Each of 

the other pairs of terms contributes a similar moving limit to the region of I 5 ,t] 
space. 

In a forthcoming work we will present a detailed discussion of the shapes and 
behavior of these correlation volumes for a series of instruments. 

5.2 Entrance slit with moving edges 

We now modify the above treatment to consider the more realistic case where the 
chopper slits move with finite speed. See O] for discussion of a related problem. 
We consider that the slit width depends on the time to, at which the beam passes 
through the slit in the following manner: 
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■ tfto< -Ti 

(,■._, a + 9toif -Ti < to < , , 

''^^">- \ a-etoifO<to<Ti ^'•^> 

ifTi<to 

where a is the maximum opening of the sht (which occurs at to = 0) and 9 = a/Ti. 
Then we rewrite equ. (p2h as 



R^n = e'^l [ " dioe^"*° r " dye'^'yy (74) 

J-T„ J-a{t„) 

= e'^// dioe'"*°^ ^ (75) 

J-T„ «Ply 

where T represents those terms in equ. ( p7[ ) which do not depend on t^ and r^. 
The last integral can be evaluated to give 



^in — 6 fc 



sin ^ (aiTp - ^i^g) sin | (aiTp + /3i^a) 
i {aiTo - PiyO) i (aiTo + PiyO) 



(76) 



The rest of the calculation proceeds as above so the only change in the result is to 
replace the sinx/a; functions depending on a^To and their associated prefactors in 
( [n] ) with the function in square brackets in (uq). The correlation function ( [Tq ) is 
seen to have the correct normalization at ai = j3iy = 0, i.e. at 5 = r = 0. The two 
terms have separated maxima occurring at 

5,±5y^-v,T^0 (77) 

where 5x is the component of 5 along the direction of the incident beam and 5y its 
component in the perpendicular direction. Thus for t = the maxima lie on the 
lines through the origin: 



-^ = ±=^-^=±=3- (78) 

Ox a Vs 

where Vs is the velocity of the slit boundary. For non-zero t the slope of the lines 
stays the same but the crossing point moves along the Sx axis with a velocity Vi . The 
correlation volume for the input beam, i. e. the overlap region where the two factors in 
( [76[ ) are significant, is thus seen to be given by a complex interplay between the slopes 
of the lines and the widths of the sin x/x functions which depend on the correlation 
lengths Syc = Li/k^a, Sxc — ti/Toki, and the correlation time r = ti/2uiiTo. 

This is a demonstration of some of the advantages of the present method, show- 
ing how a change of a component can be easily accommodated and how complex 
relationships between the energy and momentum resolution can be represented in a 
straight forward manner. 
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6 Conclusions 

Starting with Fermi's 'golden rule' for perturbation theory and assuming the inci- 
dent and final wave functions of the scattered particles to be plane waves, van Hove 
[y showed that the scattering cross section is proportional to a function of the mo- 
mentum and energy transfer, S (q, lo), which in turn is the Fourier transform of the 
time and space dependent auto-correlation function of the density fluctuations of the 
scattering system. 

However since no experiment is ever done with perfect resolution, measurements 
yield the convolution of S {q, uj) with an instrument resolution function. The deter- 
mination of the resolution and its use in the interpretation of the data are left to the 
experimenters. 

In this conventional approach a scattering experiment is divided up into two 
separate parts which are treated by two, different, mutually exclusive approximations. 
The scattering event in the sample is treated by assuming infinitely extended plane 
waves for the incident and scattered beams, while the beam is treated as consisting 
of classical point particles in order to describe its motion through the apparatus. 
This same, dual approach is also used to discuss multiple scattering: the scattering 
events are described as a scattering of plane wave states while the motion between 
two scattering points is treated as a classical trajectory. We hope to present a unified 
treatment of multiple scattering using the present approach in the near future. 

In the present work, the entire scattering experiment is described as a wave- 
optical phenomenon. This more general approach provides a rigorous justification 
for the usual dual approximation since, in most cases the correlation volumes are 
much smaller than the geometrical dimensions (beam cross sections, chopper pulse 
lengths, etc.) involved. However in some cases this is not true and the complete wave- 
optical description of the global system of instrument and scattering sample must 
be used. In addition to precision light scattering experiments such cases include 
perfect crystal neutron optics (e.g. interferometry) where the definition of the wave 
vector is so precise that the correlation lengths can be comparable to the geometrical 
dimensions. These cases are sometimes referred to as "spherical wave effects" [|9|. In 
fact these effects are natural consequences of the full wave-optical analysis presented 
here. 

Usual methods of determining the resolution include measurement of a known 
calibration sample, Monte Carlo calculations and estimates based on assuming a 
Gaussian form for the separate contributions to the overall resolution so that the 
separate contributions can be combined using the sum of squares rule. 

One method of carrying out the deconvolution is to divide the Fourier transform 
of the measured data by the Fourier transform of the resolution function. Usually 
however a model form is assumed for S{lf,oj) and some parameters of the model 
are varied until the convolution with the resolution function gives agreement with 
the measured data. 

In the present work we calculate the incoming wave function as it is shaped by 
transmission through the elements of the apparatus according to the laws of optics. 
For example on propagating through a slit and/or a chopper a wave acquires a corre- 
lation function with a finite width. After scattering, the individual waves scattered 
from a neighboring pair of points in the sample lose correlation on passing through 
the post-sample elements of the apparatus. We have shown that the probability of 
a scattered particle arriving at a detector at a given space-time location is propor- 
tional to the Fourier transform of the density fluctuation auto-correlation function 

(as shown by van Hove) but multiplied by a function HI S ,t] which represents the 
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resolution of the measurement and is the Fourier transform of the instrument resolu- 
tion as it is normally considered. This function is the product of a series of functions, 
each representing the effect of a single element of the apparatus on the measurement. 
Changes in one element can be easily accommodated by changes in the appropriate 
function. The regions of ( S ,t] space where each function is significant are the 
correlation volumes for the different elements of the apparatus. A measurement can 
be optimized by maximizing the region of overlap of all the correlation volumes. The 
fact that these volumes are seen to be moving as a function of the correlation time r, 
is a reflection of the fact that the w, and if resolutions depend on each other to some 
extent and we have an analytic description of this dependence. The present analysis 
offers a graphical method of optimizing each measurement, with the influence of all 
parameters such as the scattering angle included. In a forthcoming work we will 
give detailed examples of how this can be carried out for various types of scattering 
measurements. 

7 Appendix - Normalization and calculation of cross sections 

In the main text we have neglected all prefactors and normalizing constants in the 
Green's functions and incident wave function. In this appendix we will show how 
inclusion of these factors leads to the usual expression for the cross section and give 
a proper calculation of 4'in as mentioned in footnote 2 of section 2. We will treat the 
most general case of space and time dependent density fluctuations (sec. 5). 

7.1 Normalization of the Green's functions 
7.1.1 Calculation of the incident wave. 

To calculate the wave function incident on the scattering sample we use 

tpin{^,t)^— dto I (f'So1po{'x'o,to)lt ■^oG{~^ -l^o,t-to) (79) 

where the integral over cPSo is taken over the plane of the entrance slit, Zq = 0, 
and '0O is non-zero only in the entrance slit (—a < ijo < a). (We make the Kirchhoff 
approximation as is usual in optics) . The function G satisfies the boundary condition 
G (zo = 0) = and can be taken as 

G = .g(T^+,t)-.g(l^_,t) (80) 

where 

{^±f - (x - xl) + (y - yof + {zT Zof (81) 

and g(1^,t) satisfies 

^,^_2ma^^_^^^(3) 
in at 
g{T^,t) is then normalized as |18[ : 

._4 , /27rih\ ( m \3/2 2 
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Then (equ. ^) 



3/2 /•*+ dto 



(84) 

where G = exp ( i^r- ) is the Green's function without prefactor as used in the main 

text. 

7.1.2 Calculation of the scattered wave. 

Taking the interaction potential as 



V{T',t)^2TT—p{T',t) (85) 

m 

where b is the scattering length and p is the number density of scatterers the scattered 
wave is found to be (making the Born approximation) 

V-sc (T", t) = ^!— Go{^ -T's,t-t,)p (T's, ts) ^m {^s,ts) (fr.dts (86) 

im J 

where Go satisfies 

and is normalized as lllql 



Co-{ "^ Y^^C 



so that (86) 



where t^ ^ ^^^ the average flight times in the two arms of the apparatus and z\ is 
the distance from the entrance slit to the center of the sample. We have neglected 
the variation of the prefactors over the regions of integration as these are small 
compared to the distances and flight times involved. With suitable normalization 
of the incoming wave, '0o, the intensity, 1^, counted in the detector during the time 
interval ITd-, will be given by 



W2 / d^ri / did IV's 



V2b^ (^£j:y J^o^l^ J J d^sdts J J dHdrRout {pp') Rr. 

(90) 
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where Rout,in are defined as in section 2 and {pp') is related to the van Hove corre- 
lation function Gsi <5 , t J by [Q 

{pp')=PsGs(^,t) (91) 

and the correlation function of the input beam is defined to be 

(V.^ (T>, t) ro {^', t')).„c = So5(^) (T> - T^') 6 {t t') (92) 

I is the height of the beam perpendicular to the scattering plane. 

The multiple integral in (Q) has been evaluated in (^). Setting all the slit 
widths and chopper opening times to zero in the function H so that this takes on its 
maximum value (note: linia-^o Ji {x) jx = 1/2) we get from ( |90| ) for the intensity Id 

where we have used van Hove's definition |l|: 

S{lt,Lo)^^J J d?5dTe<^-^-^^)Gs {y, t) (94) 

N ^ Ps {jiRlt) is the total number of scattering centers in the sample and 

K = {2To) (2a) (2T,) {2Td) (2d) (95) 

Since zi = viti {— Li) we finally obtain: 



m\4 So vf 
^) {2T:fWi' 



^^-«26M^) 7ebTS3^'K5(^^-) (96) 



7.2 Calculation of the cross section 

7.2.1 Normalization of the input wave correlation function 

Rewriting the definition (|9^) we have 

= Eo^^ f d^ke'^i^-^')— I duje-'<'''') (97) 

so that 



for the white spectrum we are assuming. From this we conclude that the density 
fluctuations contained in the interval d^kduj are given by 
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diodk^ (27r) 

or 



"^ ^ d^kdcj = -^^Skduj (99) 



^'^ ^° (100) 



7.2.2 Calculation of the intensity in terms of the cross section. 

The incident flux in a velocity interval dvi can be written 

( d^p \ fdu;\ ^^^ 2 



with vi the incident velocity and dQi the input solid angle. Multiplying by the 
chopper opening time {2To) and the area of the input slit (2a/) we obtain the total 
number of particles contained in one pulse. The fraction of particles scattered into 
the velocity range dv2 and solid angle dfl2 is then given by 



where in the first factor As is the area of the sample exposed to the beam. We use 
the cross section definition of van Hove which applies to the whole sample {i.e. N 
times the cross section for the individual scattering atoms). Now with 



duj mv , vi,22Ts,d ,r^ As ,^ 2dl 

^2 



-J- - ^, dwi,2 = —^^ -, dill = -2-, dfl2 = -pr (103) 

dv h ii2 ^1 



we have 



Eq ■mviVi2Ts 2^8 ,r,rp ^ fr, ,^ f J_\ rf^cr mv2 V22Td 2dl 



for the number of scattered particles reaching tlie detector. 
Collecting terms we have 



nf \hJ till dcudn ^ ' 



{2n) 
Comparing with (Bq) we see that the two expressions agree if 



diodil vi 

which is equ. (26) of ref. |^. 



b'' — S{lt,uj) (106) 
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Figure Captions 

Fig.l) Scattering in first Born approximation (equ. pi). The incoming wave is 
scattered at T*s and the scattered wave propagates according to Go- The Green's 
function (SS) has the property that at a fixed time the wavelength decreases with 
distance from the source. 

Fig. 2 a) The incident wave field, initially uncorrelated, emerges from the entrance 
slit with correlations over a length scale d, inversely proportional to the slit width, 
a. 

2 b) Showing the phase relations responsible for the correlations of the wave 
field according to the van Cittert-Zernike theorem. Waves emerging from the uncor- 
related point sources a, b, c arrive at the scattering points r^, r^ with different relative 
phases. The phase difference of the waves arriving at the different scattering points 
from each source point depends on the location of the source point as well as the 
distance S = T^s ~ r' s, increasing with 6 . This increasing phase difference leads 
to a decrease in correlation between the waves arriving at the points ^s , r^ . Because 
the phase differences between the uncorrelated source points are random they have 
no influence on the correlation function. 

Fig. 3) Illustrating equ. (|l^). The functions Rm, Rout represent the effects of the 
phase differences between the paths 1 and 2 on the correlation functions of the beam. 
Gs is the correlation of the density fluctuations at the different scattering points and 
l^dl is the sum of the contributions from paths going through all possible pairs of 
points in the sample. Only pairs whose separation lies within the correlation volumes 
contribute. 

Fig. 4) Elastic scattering - general case. Showing the meaning of the symbols used 
in the text. S is the distance between the two scattering events at T*s, T*^. 'ts, "e*s 
are the positions relative to center of the sample. Further symbols are defined in the 
text. 

Fig. 5) Resolution function for small angle scattering 

Fig. 6) Scattering from a time dependent system. The scatterer is assumed to 
be concentrated at a single point (~r*s)- The resolution is determined by the phase 
difference between paths (1) and (2). Note the analogy with elastic scattering - fig. 
4. T is the delay between the two scattering events at ts,t'g. S,s,£.s ^^^ the delays 
relative to the nominal scattering time is- 
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